Production of Gravitational Waves in the nMSSM 

Stephan J. Huber*^^^ and Thomas Konstandin'^^l*] 

(1) Department of Physics and Astronomy, 
University of Sussex, Brighton, BNl 9QJ, UK and 
(2) Department of Theoretical Physics, 

OO ' 

' Royal Institute of Technology (KTH), AlbaNova University Center, 

■ Roslagstullsbacken 21, 106 91 Stockholm, Sweden 

(Dated: May 28, 2008) 

^: 

During a strongly first-order phase transition gravitational waves are produced by 

OO . 

CN . bubble collisions and turbulent plasma motion. We analyze the relevant character- 



X 



istics of the electroweak phase transition in the nMSSM to determine the generated 
gravitational wave signal. Additionally, we comment on correlations between the 



^ . production of gravitational waves and baryogenesis. We conclude that the gravi- 

tational wave relic density in this model is generically too small to be detected in 

(N 

J> ' the near future by the LISA experiment. We also consider the case of a "Standard 

, of gravitational waves. 

o 
o 



I. INTRODUCTION 

Presently, several experiments are under consideration that could detect for the first time 
a stochastic background of gravitational waves (GWs). One of the main motivations to 
pursue these experiments is that the discovery of a relic gravitational background would be 
a smoking gun signal from inflation [l| and hence might allow to test the paradigm of an 
era of exponential expansion in the early Universe. 

Another source of stochastic GWs are strongly first-order phase transitions. Space-based 
experiments, such as LISA [2] and BBO [3] will have remarkably good sensitivity at fre- 
quencies that coincide with the redshifted spectrum of GWs produced during an electroweak 
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phase transition at temperatures T ~ 100 GeV. This opens the possibihty to infer informa- 
tion about the electroweak phase transition from GW observations. A strongly first-order 
phase transition not only would produce GWs, but could also generate the observed baryon 
asymmetry of the Universe (BAU) through electroweak baryogenesis 4], [Sj]. This way the ob- 
servation of a stochastic background of GWs might teach us about electroweak baryogenesis 
in extensions of the Standard Model (SM). 

During a first-order phase transition there are two distinct mechanisms that produce 
GWs: The colliding phase boundaries [gI, 3, Isl, and the turbulent motion of the plasma 
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10, 
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Il3l . Il4| . Concerning the latter, there is an ongoing discussion in the literature 
on how to correctly model the temporal correlation of the turbulent plasma. In particular, 
the different approaches predict different peak frequencies of the resulting GW spectrum. 
Studies of GWs from the electroweak phase transition have been performed in Refs. hi 

n 

and pJ!] without relying on particular particle physics models. The results are presented 
as functions of the two main parameters of the problem, which are the typical size of the 
colliding bubbles and the available energy. It was concluded that a sufficiently strong phase 
transition could lead to an observable GW signal at LISA. Of course, in a given model the 
bubble size and the available energy are linked, and it is interesting so see what happens to 
the GW signal. The cases of the MSSM and the NMSSM (Next-to-minimal supersymmetric 
SM) were studied in Ref. [l3]. While it was found that in the MSSM the produced amount 
GWs is orders of magnitude below the LISA sensitivity, the situation in the NMSSM seems 
much more promising. The cubic terms present in the NMSSM tree-level Higgs potential 
can lead to a much stronger phase transition 18|. 



However, the authors of Ref. 



17| used a crude method to determine the bubble configura- 



tions, which overestimates the strength of the phase transition and the GW signal. The aim 
of this work is to improve on this point to arrive at a more realistic estimate of the signal 
strength. To achieve this we apply our recently presented numerical method to compute the 
bubble shapes jl9| . 



Like in Ref. [l7j we consider a singlet extension of the MSSM, in our case the nearly 
Minimal Su persy mmetric Standard Model (nMSSM) 2c]]. In contrast to the models studied 



in Refs. 



17, 



18| . it solves the yU-problem of the MSSM, without creating a domain wall 



problem or destabilizing the electroweak hierarchy through divergences of the singlet tadpole. 
With this motivation, the cosmology of the nMSSM has been been studied intensively in 
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the literature. It was shown in Ref. 



2l| that the model can have a strong enough phase 



transition for baryogenesis, while the lightest neutralino can at the same time provide the 
dark matter of the Universe. The baryon production was studied with a positive result 
recently in Ref. [2^. Finally, Ref. 23| investigated to what extent these cosmological issues 
can be studied at colliders, in particular at the LHC 

In this paper, we compute the GW signal in the nMSSM using numerically determined 
bubble configurations. In addition, correlations between the amount of the produced GWs 
and the generated BAU are discussed. Our results for GW production are less optimistic 
than the estimates of Ref. 17|], making it doubtful for LISA to detect a GW signal, even for 
an optimal choice of parameters. BBO, on the other hand, could detect GWs if the nMSSM 
phase transition is extremely strong. This study is a continuation of work on the baryon 



asymmetry in the nMSSM that has been published recently in Ref. 



22] . Since the formalism 



to determine the BAU, the analysis of the phase transition and phenomenological aspects of 



the nMSSM are discussed in detail in Refs. 



21 



2211 ■ we will be rather brief on these issues. We 



also study GW production in the SM augmented by dimension-six operators. Such a model 
can be viewed as the effective low energy description of some strongly coupled dynamics at 
the TeV scale. As this model has only one Higgs field and few parameters, the analysis is 
much simplified and more transparent compared to the nMSSM case. It has been shown 



24 



26|. 



that an if^-term in the Higgs potential easily leads to a strong phase transition 
Additional dimension-six operators can also provide new sources of CP violation that can 
account for the observed baryon asymmetry {3, Q without generating large electric dipole 
moments 



The paper is organized as follows. The next section contains the formalism to determine 
the relevant properties of the phase transition. In Sec. IIIII the model is discussed, while 
in Sec. llVl the nMSSM and its temperature dependent effective potential is summarized. In 
Sec. |V] the numerical results for the nMSSM are presented before we conclude in Sec. IVIII 
Appendix El contains a brief review of first-order phase transitions in cosmology. 



II. GRAVITATIONAL WAVES 



The discussion of GW production in this section closely follows the analysis in Ref. 15|. 
Details about the GW production during a first-order phase transition can be found in this 



4 



work and the references therein. 

One source of GWs during the phase transition resuhs from coUiding Higgs bubbles at 
the end of the phase transition. Depending on the strength of the phase transition, the 
bubble wall profile propagates faster or slower than the speed of sound, which is l/-\/3 in 
a relativistic thermal bath. In the former case GWs are produced by detonation, while in 
the latter case deflagration is the dominating process. Since GW production is strongly 
suppressed in a, deflagration, we focus in the following on the detonation mode as discussed 
in Refs. {g, ^ 



m 



29|. 



To determine the magnitude of the produced GWs, two parameters of the phase transition 
are essential. The first one is related to the latent heat, 



dV 



(1) 

which provides the energy available to be transferred to GWs. The amount of GWs depends 
on the ratio of the latent heat to the energy density of the radiation in the plasma, 

This parameter measures the strength of the phase transition, as relevant for the GW pro- 
duction. T^, denotes the temperature of the phase transition, and g^, is the effective number 
of degrees of freedom (in the SM: g^ = 107.75). There are several prescriptions to define 
the start and end of the phase transition, and the time of maximal GW production. This 
temperature is further specified in the Appendix, see Eq. ( lAlOp . 

The second parameter that enters the production of GWs is the typical radius of the 
colliding bubbles, (R). It sets the length scale of the problem. It can roughly be expressed 
as the product of the wall velocity, v^, and the duration of the phase transition, r. The 
latter can be approximately related to the logarithmic time derivative, /3, of the nucleation 
rate, r ~ f3~^. In adiabatic approximation one has 



dT V T 



(3) 



where H^, denotes the Hubble parameter at T*. Hence, the typical length scale is proportional 
to 

{R)(xVbT^j. (4) 



In contrast to Ref . [17[ , we state our results in terms of the typical radius of the detonating 
bubbles (R) instead of using the approximate expression containing jS. It is not known what 



is the optimal definition for (R). Following Ref. 19], we use the maximum of the bubble 
volume distribution, which should be closely related to the typical scale for energy injection 
into the turbulent plasma. In this case, one obtains 



(5) 



The relevant formulas can be found in the Appendix. 

The energy density in GWs per logarithmic unit of frequency is usually normalized to 
the critical density, 

0(/) = _L%^. (6) 



Pcrit d\nf 

The GW contribution from colliding bubbles at the peak frequency, as observed today, is 
then given by ol 



1.2 X 10- V {RfH^, 



a 



1 2 



a + 1 



Vb 



0.24 + 



100 



nl/3 



(7) 



Assuming a detonation, i.e. a strong phase transition with a supersonic bubble wall, the wall 
velocity Vb is approximately given by 291] 

1/ V3 + \/«^ + 2a/3 



Vb{a) 



1 + a 



(8) 



This is a good approximation for very strong phase transitions. Finally, the efficiency 



K is 



1 



1 + 0.715a 



0.715a + 




4 /3a 
27 



(9) 



It measures the fraction of latent heat that is transformed into bulk motion of the fluid. 

When the 



The second source of GWs during the phase transition is turbulence j^, Q, Q • 
bubbles collide, the plasma is stirred up, which creates a cascade of eddies in the plasma. The 
relic GW density at the peak frequency coming from turbulence, according to Refs. |l^, 
is found to be 



15|, 



"'0'^ 'turb 



5.6 X IO-'^mI {RfH^ 



100 

g* 



1/3 



(10) 



As additional parameter, the characteristic turnover velocity us of the eddies in the plasma 
enters, given by 15] 
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The expression ( ITO|) relies on a certain modeling of time correlations in the turbulence, which 



was introduced in Ref. 



111 ]. Similar results were obtained using Kraichnan's ansatz for the 



time correlations 14 1. 

Recently, a different approach in calculating the contribution from turbulence was pre- 
sented in Ref. ISj- This work is based on Richardson's model of turbulent motion which 



13| emphasize 



is based on random velocities of the turbulent fluid. The authors of Ref. 
the fact that the GWs inherit the momentum spectrum from the eddies rather than their 
frequency spectrum. The resulting GW density at the peak frequency is 

M 1/3 



/io^Richardson ^ 6.7 X IQ-^^g {Ry H, 



100 

9* 



(12) 



Notice that if the results are stated in terms of (R), as done here, the numerical coefficients 
appearing in the GW densities (and frequencies) of the two approaches to turbulence ba- 
sically agree. As mentioned before and detailed in Ref. the main difference is based 
on the question if the GWs inherit the momentum spectrum or the frequency spectrum of 
the eddies resulting into an addition factor us in the GW density (and also in the peak 
frequency) and different scaling behaviors. 

On the other hand, comparing the numerical coefficients for turbulence in terms of (3 as it 
was presented in Ref. 12] and in the appendix of Ref. ISj one observes a disagreement that 
can be traced back to the fact that different expressions for relating the typical length scale 
{R) to P have beens used. In Ref. [l3| the typical bubble radius is given by {R) ~ Vb/P, while 



in Ref. [l2j the relation (R) ^ 5vb/f3 was used. To relate (3 to (R) we use for both approaches 
the relation given in Eq. what is motivated in the appendix. Deciding between these 
different approaches to turbulence is beyond the scope of the current paper. We simply 
present results for both cases, indicating the spread in the theoretical predictions. 

For weaker phase transitions the contributions from bubble collisions are larger than the 
ones from turbulence, but if the phase transition is really strong, us ^ Vh, both contributions 
can be of comparable magnitude or the latter can even dominate. 

Since we assume the phase transition to lead to supersonic bubble wall expansion and 
use the GW production formula at the frequency peak, the given production rates of GWs 
should be understood as upper bounds. If the wall velocity is subsonic, detonation does 
not take place and the contribution of turbulence is strongly reduced, not only because the 
velocity Vh is small, but also because the efficiency factor k, is much smaller than in Eq. IQ, 
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which as well assumes supersonic bubble velocities. 

To roughly determine the spectrum of the GWs, it is sufficient to know the peak frequency 
and the scaling behavior. The peak frequencies are given by [9 

1 n 



12 



13| 



dct 



fi 



turb 



/] 



Richardson 



1.6 X IQ-^mRzVb 

1.7 X IQ-^mRzus 
1.7 X lO^^mHz 



(R) 100 GeV 
1 



9* 



1/6 



{R)H, 100 GeV 



T 



9^ 



100 

" 9* 
.lOOJ 

1/6 



1/6 



LlOOJ 



(13) 
(14) 
(15) 



{R) 100 GeV 

Notice that the peak frequency in Richardson's model differs from the "usual" result for a 
turbulent plasma by a factor us due to the dispersion relation of the eddies. This shifts the 
peak frequency to higher values. This fact is relevant observationally, as the peak frequencies 
for GWs from the electroweak phase transition are typically at the lower end of the LISA 
window, where its sensitivity is already reduced. 

The GW spectrum from bubble collisions rises as /^'^ below and decreases as f~^'^ above 
the peak [g^. The corresponding scaling behavior in the case of turbulence is given by /^ ^ and 



/ respectively 
model is given by 



Finally, the scaling behavior in the analysis based on Richardson's 



^0 ^Richardson OC < 



/^ / < 2M5/peak 

/; 2^5 /peak < / < /peak 
/ ) /peak ^ / 



(16) 



III. THE STANDARD MODEL WITH DIMENSION-SIX OPERATORS 



Before turning to the nMSSM, let us discuss GW production in a simple "toy model". 
It consists of the SM augmented by certain dimension-six operators. These operators 
parametrize unknown physics at the cutoff scale, typically at a few hundred GeV. They 
could originate from some new strong dynamics, such as technicolor or weak scale gravity, 
or from simply integrating out scalar fields. 

The model contains a single Higgs doublet, H, whose potential is stabilized by an 



interaction 



24 



26| 



rr2 



X 



^ ^ 2 4 



1 



(17) 
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This potential has two free parameters, the suppression scale of the dimension-six operator, 
M, and the quartic coupling, A. The latter can be eliminated in terms of the physical Higgs 
mass, rriH- Since the potential is stabilized by the term, A can be negative. In this 
case a barrier in the Higgs potential is present at tree-level, which induces a strong first- 
order electroweak phase transition. Evaluating the one-loop thermal potential, it was shown 
m Ref. 126^ that thephase transition is strong enough to avoid baryon number washout, 
i.e. {H)tJTc > 1.1 if M < 850 GeV and = 115 GeV. At the critical temperature, 
Tc, the broken and symmetric phases are degenerate in energy. Taking M = 500 GeV, a 
strong phase transition is present for itlh < 180 GeV. 

Dimension-six operators also induce new sources of CP violation. In addition 
to the ordinary Yukawa interaction of the top quark, ytHf^qs, one has an operator 
{xt/ M'^){W H) Hf^q^. Along the bubble wall these two contributions to the top quark mass 
enter with varying weight. In this way a possible relative phase between the two operators 
induces a varying complex phase in the top quark mass such that tops and anti-tops behave 
differently in the bubble background. Chiral charges are built up in front of the bubble wall, 
which the sphalerons transform into a baryon asymmetry. A semi-classical analysis of the 



these processes has been performed in Refs. 



26 



271]. It was shown that the observed BAU 



can indeed be produced this way. This can be achieved without generating dangerously 



large electric dipole moments 



28|. 



To compute the spectrum of GWs one has to study the real time history of the electroweak 
phase transition, in particular one must determine the parameters a and {R). This requires 
knowledge of the bubble nucleation rate, which depends on the energy of the nucleating 
bubbles and varies with temperature. The bubble configuration is a stationary solution to 
the Euclidean equations of motion of the Higgs fields. The relevant formulas are collected 
in the Appendix. We use the one-loop thermal Higgs potential as described in Ref. {2^] but 
evaluate the full one-loop result instead of using the high-temperature expansion. Since there 
is only one Higgs field, the bubble solution can be obtained by a simple over/under-shooting 
method. 

We define the onset of the phase transition as the temperature, T„, where the first bubble 
per horizon volume nucleates, (lA7p . As the Universe cools down more bubbles nucleate, 
expand and finally collide. During this stage one can neglect the expansion of the Universe, 
which should only matter for extremely strong phase transitions, /3 ~ iJ*. Most bubbles 
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FIG. 1: Different characteristics of the phase transition as functions of the parameter M. The 
Higgs mass is chosen to be m/j = 120 GeV. 



collide at the end of the phase transition, which can be defined as the temperature, Tf where 
a fraction of 1 — 1/e = 0.632 of the volume is converted into the broken phase, (lAlOp . Finally, 
we obtain (R) from the maximum of the bubble volume distribution at Tf, as mentioned 
before. 

In our numerical example we use the Higgs mass mn = 120. We have checked that 
for larger Higgs masses the results are essentially the same, only the relevant values of M 
change. When M decreases, the the strength of the phase transition increases. At M ~ 576 
GeV the system becomes metastable and the Universe gets stuck in the symmetric phase. 
As shown in Fig. [T] the phase transition is strong enough to avoid washout after electroweak 
baryogenesis. While Tc and (H) /T^ change only slowly in the presented range of M, Tf and 
T„ drop rapidly. This signals that the system is close to metastability. At the same time 
a rises from a value around 0.1 to order unity and the typical bubble radius {R) increases 
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f/Hz f/Hz 
FIG. 2: The spectrum of the density of GWs for different values of the parameter M. The used 
parameters are given in Tab. HI Smaller values of M lead hereby to a stronger phase transition 
and larger GW production. In the shaded region, the sensitivity of LISA and BBO is expected to 
drop considerably. In the left (right) panel, the turbulent contribution to the GW spectrum from 
Ref. Q (Q) has been used. 



set 


M 1 GeV 


a 




T, / GeV 


1 


600.0 


0.128 


0.015 


62.8 


2 


588.0 


0.201 


0.028 


53.1 


3 


582.0 


0.311 


0.050 


45.6 


4 


578.0 


0.586 


0.116 


37.4 


5 


576.5 


1.197 


0.318 


30.4 


6 


576.0 


2.268 


0.592 


25.6 



TABLE I: Sets of parameters used in Fig. [2j 



dramatically. In Fig. [2] the GW spectrum for several choices of M is shown, as given in 
Tab. [H together with the LISA and BBO sensitivities. In gray the region of / < 10~^ Hz 
is indicated, where observational sensitivities drop considerably. In the two panels, the two 
different predictions for GW production from turbulence as given in Refs. 12] and [3] have 
been used. The plot demonstrates that as the phase transition becomes stronger and the 
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amplitude of the GW signal becomes larger, the peak frequency moves to lower values. This 
is directly related to the smaller values Tf and larger values of {R)H^: in Fig. [H As a result, 
even in the case of an extremely strong transition, M ~ 576 GeV, it will be difficult for 
LISA to detect a GW signal. This would require sensitivity at lower frequencies. Things are 
different for BBO. It would be able to detect GWs if M < 585 GeV, which corresponds to a 
few percent tuning in M. In this case also the subtraction of backgrounds, in particular from 
white dwarf binaries, has to be under control (see |l^, and references therein). The main 
difference between the two models of turbulence is that in the case of the Richardson model 
the two-peak structure in the GW spectrum disappears, as the turbulence contribution is 
shifted to higher frequencies. If present and observed, the double-peak structure could help 
to identify a phase transition as source of the GW signal. Comparing with the general 
analysis in Ref. [16|] we arrive at a more pessimistic picture because large values of a mean 
small (R) and a low peak frequency. The decrease in temperature for very strong phase 
transitions additionally enhances this effect. 



IV. THE NMSSM 



While the electroweak phase transition in the MSSM is too weak to produce observ- 
able GWs 171], the strength of the phase transition increases considerably if the MSSM is 
equipped with an add.tiona. gauge singlet Q. There are various possibihti^ to .mple- 
ment this idea. A phenomenologically attractive candidate is the nMSSM [20| with the 
superpotential 



W^nMSSM — XSHi ■ H2 



m 



12 



A 



S + Wi 



MSSM- 



(18) 



In contrast to the models studied in Ref. 18|, it solves the /x-problem of the MSSM without 
destabilizing the electroweak hierarchy or the generation of domain walls at the electroweak 
phase transition. However, the discrete R-symmetry necessary to accomplish this task, for- 
bids the singlet self-couphng and leads to a quite constrained Higgs and neutralino/chargino 



phenomenology 



21 



22|. In particular, the nMSSM does not contain a cubic singlet term in 



the superpotential, unlike the NMSSM. Hence, the singlino obtains its mass only by mixing 
with the Higgsinos. It is typically very light, and can be a dark matter candidate 211]. An 
early version of a singlet model without a self-coupling can be found in Ref. 32]. 
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The tree-level Higgs potential of the nMSSM including soft SUSY breaking terms reads 



V = Vf + Vd + V^oft, 



where 



Vf 
Vd 



X'\Hi ■ H2\'' + \^\S\\hIHi + HIH2) - {mi^Hi ■ H2 + h.c), 



9' + 9" 



{HIH2 - HlHif + ^\HIH2 



mlHlHi + mlHlH2 + ml\S\^ + (tsS + h.c.) + {axSHi ■ H2 + h.c). 



(19) 

(20) 
(21) 

(22) 



Additionally, we take into account the Coleman- Weinberg one-loop terms 

1 



167r2 



^gbh{ml) -^g}h{m)) 



(23) 



b f 

where the two sums run over bosons and fermions with the degrees of freedom Qf, and Qf 
respectively, and 

,4 



h[m 



m 
T 





3" 


MS) 






~ 2_ 



(24) 



At finite temperature, the effective potential is modified by the following thermal one-loop 
contributions: 

rp4 



27r2 



Y,9Mrnl/T') -J29fJ-H/T') 

■ b f 



(25) 



with 



/•oo 

My^)= dxx^ log(lTexp(-v/x2 + y2)). (26) 
Jo 

For both one-loop contributions we use the following numbers of degrees of freedom 

gw = 6, gz = 3, gt = 12, g^^ = gi^ = 6, (27) 



and the renormalization scale is chosen to be Q = 150 GeV. As in Ref . 22] , the stop masses 
are chosen as 

mu = mQ = 500 GeV, (28) 

such that the strength of the phase transition originates mostly from the (tree level) singlet 
sector. As most SUSY particles are assumed to be heavy, we approximate the number of 
effective degrees of freedom by the SM value, g^ = 107.75. 
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At zero temperature, both Higgs fields and the singlet field acquire a vev and the Higgs 
vevs are constrained by 

0(T = 0) = t; ~ 174 GeV, (l)\T) = \ p + | (H^,)^ \\ (29) 

Additionally, since one of the chargino masses is approximately given by 

m^± ^ A(^) > 114GeV, (30) 

the parameter space is restricted to regions with a sufficiently large singlet vev. The lower 
bounds on the masses of the physical Higgs states further restrict the parameter space, 
such that the nMSSM constitutes a much stronger constrained model than models with an 
explicit Higgsino/chargino mass (/i-term) or a singlet self-coupling. These bounds are easier 
to satisfy for larger values of A, which therefore prefers to be close (even beyond) the Landau 
pole value. 



V. NUMERICAL ANALYSIS 



Our numerical analysis follows the approach already used in Ref. 22] to determine the 
BAU in the nMSSM. The relevant eight free parameters of the nMSSM are chosen randomly. 
The generated parameter sets are then confronted with several constraints on the mass 



spectrum and the Z-width (for details see Ref. [22] )• In addition, we restrict ourselves to 
cases which promise a rather strong phase transition, 0(Tc)/Tc > 1, by inspection of the 
potential. If the parameter set passes these constraints, the properties of the phase transition 
and the GW relic density are determined. 

□ 

Compared to the former work on the GW production in the NMSSM [17], the analysis is 
improved in the following points. First, we use up-to-date bounds on the particle spectrum, 
which severely constrains the parameter space. Second, the bubble configurations of the six 



scalar fields are determined exactly, while in Ref. [17| this problem was reduced effectively to 
a one-dimensional one, which overestimates the strength of the phase transition and hence 
the GW production. Third, the parameters relevant for baryogenesis are determined in the 
specific cases and hence a direct correlation between the BAU and the magnitude of the 
GWs can be inferred. 

Because of the linear and tri-linear terms of the singlet in the potential (l20i) . a stro ngly 



first-order phase transition in the nMSSM is possible due to tree-level dynamics 
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FIG. 3: The left panel shows the correlation between the parameters a and the typical size of the 
bubbles at the end of the phase transition, (R). The right panel compares the typical bubble size 



with P at first nucleation, T = Tn, and the end of the phase transition, T = T 



34( 1 ■ The phase transition in this case is described by six fields: The two vevs of the neutral 



Higgs fields {H^), {H2), the vev of the singlet (S) and their complex phases; however, the 
relative phase between the Higgs vevs does not enter the potential and can be removed using 
its equation of motion. In this work, the profiles of the vevs during the phase transition are 
calculated numerically using the method presented in Ref. R that is briefly described in 
Appendix \M 

The left panel of fig. [3] shows the correlation between the parameters a, the available 
energy, and (R) for approximately 150 parameter sets. We find the same behavior as 
in the SM with dimension-six operators. The right panel compares the typical bubble size 
with (3 at first nucleation, T = Tn-, and the end of the phase transition, T = Tf. Note that 
/3(T„) usually underestimates the bubble radius, while P(Tf) overestimates it. For strong 
phase transitions, f3 varies considerably between T„ and Tf. This makes our more careful 
determination of (R) necessary. The same applies to the SM with dimension-six operators. 

Since there is a relatively strong correlation between a and the bubble radius, we correlate 
in the following only the parameter a with the parameters that enter the determination of 
the BAU and the production of GWs. 

Figure H] shows several characteristic properties of the phase transition. In all parameter 
sets under consideration, baryon number washout after the phase transition is sufficiently 
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FIG. 4: The correlation between the parameter a and several properties of the phase transition: 
The temperature at first nucleation T„, the temperature at the end of the phase transition Tj, the 
thickness of the bubble wall 1^ and the ratio between the Higgs vev and the temperature 4>{Tf) /Tj. 



suppressed [31|, 4>(Tf)/Tf > 1.1. As concluded in Ref. [221], electroweak baryogenesis is a 
generic feature for this class of models. However, this statement is based on the assumption 
that the wall velocity is subsonic, such that diffusion is operative. From the parameters that 
enter the generation of the BAU in Eq. (l38ll . three are expected to have a strong dependence 
on the strength of the phase transition: The wall thickness and the temperatures Tn and 
Tf. These quantities are also plotted against the parameter a in Fig. HI Both, the wall 
thickness and the temperatures, decrease with increasing strength of the phase transition. 

The magnitude of the produced gravitational waves is displayed in Fig. [5], separately for 
the contribution from collisions and turbulence. For our parameter sets, the strength of the 
phase transition never exceeds a = 0.2 and the density of the produced GWs is too low to 
be observed by the LISA experiment that has a sensitivity of HqQ ~ 10"^^. On the other 
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hand, the BBO experiment has a maximal sensitivity of h^Q ~ 10^^'' which is several orders 
below the maximal GW density we found, h^Q ~ 10^^^. 

However, even the prospects for the BBO experiment to observe GWs from a strongly 
first-order phase transition in the nMSSM are very limited as seen in Fig. [61 The plot 
shows six spectra corresponding to the parameters given in Table [Tll Even for relatively 
strong phase transitions, the spectrum of GWs would escape detection, since a stronger 
phase transition proceeds at lower temperatures and leads to larger bubbles, as can be seen 
in Fig. m This effect lowers the peak frequencies according to Eqs. (|T3l) to regions where the 
sensitivity of BBO is strongly reduced. This is the same phenomenon we already observed 
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set 
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T, / GeV 
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0.03 


0.003 


130 
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0.05 


0.01 


110 
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0.07 


0.03 


85 
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0.1 


0.1 


80 
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0.15 


0.1 


70 


6 


0.2 


0.1 


60 



TABLE IL Sets of parameters used in Fig. [H 



in the SM with dimension-six operators. 

This leaves the question why the nMSSM does not seem to have a region in parameter 
space where a is very large, as observed in the NMSSM in Ref. 17|. One might think 
that the nMSSM is much more constrained than the NMSSM and that Higgs and chargino 
masses in accordance with current LEP bounds strongly restrict the parameter space of 
viable models. Nevertheless, the nMSSM admits very strong phase transitions induced by 
tree-level terms, even leading to met ast ability. So the crucial difference seems to be that we 
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determine the tunnel action exactly, while in the work [17|] the problem was reduced to one 
dimension, which can considerably overestimate the strength of the phase transition. 

According to Eq. ([2]), large a requires either large amounts of latent heat or small nucle- 
ation temperatures. The latent heat cannot be made arbitrarily large. For example, very 
small/large values for a\ or ts can lead to an unbounded effective one- loop potential or 



are not compatible with the mass constraints (one explicit example is given in Ref. |22| ) . 
Thus, the upper bound of the latent heat is expected to be limited approximately by the 
electroweak scale (f ~ 1 74 GeV) and stronger phase transitions require smaller temperatures 

o^<M^t^• (31) 



On the other hand, the nucleation temperature T* cannot be many orders smaller than the 
temperature T^, where the two minima of the effective potential are degenerate. This can 
be seen in the following way. Assume that the Higgs vev is close to its zero temperature 
value, 0(T) ti, and that the potential is of the form 

A\/(T) = Al^(O) + f = f (T^ - T^) v^, (32) 

with some constant g that mostly depends on the particle content of the theory. Since the 
three-dimensional action 5*3 scales as AV""^ for small AV, the condition in Eq. flA13|) turns 
into 



T 



T.' (T,^-T.T ^''^ 

This implies > |Tc, since the critical temperature has to be larger than the minimum 
of this function. The numerical results in the nMSSM show that the minimum is rather 
given by T^, > cTc with the constant c ~ 0.5. This is mainly due to the fact that the 
three-dimensional action only scales as AV'"^ for very small AV. The closer T* is to the 
minimum, the larger is a and the smaller is (3. This in turn leads to the fact that small 
critical temperatures require small latent heat, since one obtains, using Eq. fl32l) . 

\AV{0) \<e, = f (T^ + T^) v'<{l + c~^) f v\ (34) 

Thus, the parameter a scales in this limit as 



Ti AV T3 
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Moreover, small latent heat e^, and thus small |AI^(0)| tend to increase the three- 
dimensional tunnel action 

140 ^ S^{T,) oc AV^(T,)-2 > AV^(O)-^ (36) 

such that the condition in Eq. flAlSp implies that large a not only requires a small |Ay(0)|, 
but also a very small potential well between the two minima. Our numerical analysis indi- 
cates that this seems not to be compatible with the mass constraints on the charginos and 
Higgs bosons in the nMSSM, what limits the strength of the phase transition. Nevertheless, 
with further tuning a somewhat larger value for a might be possible. 



VI. ELECTROWEAK BARYOGENESIS 



While the plasma is driven out of equilibrium by the expansion of the Higgs bubbles during 
the phase transition, CP violation is provided by the complex coefficients in the potential of 
Eq. fl20|) . After redefinition of the fields, all sources of CP violation can be attributed solely 
to the parameter tg by a complex phase. (We do not consider a CP-violating phase in the 
Wino mass, as was done in Ref. {21].) As mentioned before, the nMSSM does not contain an 
explicit Higgsino mass /i, but only an induced Higgsino mass from the singlet-Higgs coupling 
according to 

/i = A (^) . (37) 

The change in the complex phase of the singlet during the phase transition gives rise to 
a CP-violating current of charginos in the bubble wall. Via Yukawa interactions, this CP 
violation is communicated to the tops that bias the sphaleron process. This mechanism 
is called chargino mediated electroweak baryogenesis and was analyzed in the MSSM in 



Refs. 
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36 



37, 



series of papers 3^ ^'^^ recently applied to the nMSSM 

In the case of a large Wino mass parameter, M2 > 1 TeV, the change in the complex 
phase, AArg(/i), during the phase transition is the main source of baryon number generation 
and a good estimate of the predicted 7710 is given by 

3 

r/10 ^ c{T^) ^— ^ Hf ~^ exp(-/io/rT,), (38) 



. The formalism we use to determine the BAU was developed in the 
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FIG. 7: The baryon asymmetry 7710 for different critical temperatures T^, AArg(;u) 
10. 



vr/lO and 



/ T 



where /iq denotes the parameter in the symmetric phase, A/x the change in during 
the phase transition and 1^ the thickness of the bubble wall. The two coefficients are 
c(T*) ^ 1.6T*/GeV and r ^ 0.78. The function r^io is plotted in Fig. [7] for /i = A/i and 
several values of and shows very good agreement with the full numerical results of the 
diffusion equations given in Ref. 22j. This has to be compared with the observed baryon 



asymmetry that, normalized to the entropy density, is given by 

riB - riB 
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■10 



10 



0.87 ±0.03. 



(39) 



We see that under these assumptions, electroweak baryogenesis is a quite generic feature of 
the nMSSM. However, these results are not directly applicable to the present case, since they 
have been obtained using a wall velocity smaller than the speed of sound. In this regime, 
the produced BAU does not strongly depend on the precise velocity of the wall. In the 



context of baryogenesis and for not too strong phase transitions, subsonic bubble ve 
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ocities 



are motivated by calculations of the friction effects acting on the bubble wall 

For the case at hand, the phase transition is usually strong enough to make friction 
negligible such that the phase transition indeed proceeds by detonation and with supersonic 



wall velocities, as can be estimated in the following way. According to Ref. 



4l[ |. the main 



21 



friction on the wall constitutes the interaction of the W bosons and is given by 

^ ~ log(mH-L), (40) 

where mw and rriD denote the mass and Debye mass of the W bosons in the plasma. On 
the other hand, the expansion of the bubbles is driven by the latent heat. The resulting 
pressure difference acting on the wall can be estimated to be 

Thus, in the present case friction from the W bosons can be neglected, r] ^ Ap, as long as 
a ^ 0.01 and the phase transition proceeds by detonation. 

As has been explained earlier, for an observable production of GWs, wall velocities beyond 
the speed of sound are necessary. This is problematic for standard electroweak baryogenesis, 
since in this regime diffusion processes are strongly suppressed. However, even if the analysis 
based on diffusion is not valid, other mechanisms, as e.g. quantum mechanical reflection at 
the wall |43!] might be responsible for the production of the baryon asymmetry. The relevant 
parameters that enter a quantitative analysis of the produced BAU are in fact similar. For 
example, in the case of the nMSSM with a very strong phase transition, one would expect 
that the produced baryon asymmetry is small since the low temperature leads to Boltzmann 
suppressed chargino densities. This problem might be avoided in scenarios where the CP 
violation is supplied by SM particles that are massless in the symmetric phase, as e.g. in the 



271] . Nevertheless, a strongly first- 



model with dimension-six operators discussed in Ref. [26|, 
order phase transition with supersonic bubble velocities seems to disfavor baryon production 
at the electroweak scale. 



VII. CONCLUSIONS 



In conclusion, we found that the density of GWs produced during a strongly first-order 
electroweak phase transition is in a large portion of parameter space too small to be de- 
tectable by LISA or even the BBO experiment, both for the nMSSM and the Standard 
Model with dimension-six operators. This is partially because stronger phase transitions 
proceed at lower temperatures and create larger bubbles, which shifts the peak of the spec- 
trum of the GW density to lower frequencies. This correlation was not taken into account 
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in the model independent analysis of Ref . 16| , which therefore led to more optimistic con- 
clusions. For the nMSSM our results are best summarized by Fig. [6] which displays the 
spectrum of GW densities for several parameter sets, and for the Standard Model with 
dimension-six operators by Fig. [21 Moreover, we observe that very strong phase transitions, 
a > 0.2, in the nMSSM are not compatible with constraints on the particle spectrum, or 
at least require a significant amount of tuning. The situation in the Standard Model with 
dimension-six operators is slightly more promising. Besides, it is easier to understand the 
effects of parameter tuning in this case. Generating observable GWs at BBO requires a few 
percent tuning in the suppression scale of the dimension-six operator. 

Technically, our main improvement was the full numerical computation of the bubble 
configurations, which leads to more realistic predictions on the properties of the phase 
transition that enter the GW computation. The simple approximation to bubble profiles in 
Ref. [17| seems to considerably overestimate the strength of the phase transition and in turn 
the GW signal. 

On electroweak scales, a strongly first-order phase transition is a necessary prerequisite 
for electroweak baryogenesis as well as for sizable GW production. However, there is no pos- 
itive correlation between GW observations and electroweak baryogenesis. In the nMSSM, 
for all parameter sets under consideration, washout is suppressed after the phase transition 
and electroweak baryogenesis might be possible. Besides, GW production is marginal for 
many parameter sets, such that electroweak baryogenesis cannot be excluded, even if GWs 
are not observed down to a sensitivity of h^Q = 10^^'^. On the other hand, observing GWs 
in the near future does not directly support electroweak baryogenesis. In fact, the oppo- 
site might be true. Since significant GW production by collisions and turbulence requires 
supersonic bubble wall expansion, diffusion is suppressed which reduces the produced BAU 
in mechanisms based on transport. Due to rather thin walls, other mechanisms, as e.g. 
CP-violating reflection, might be effective in this regime, but reliable predictions in this 
scenario are missing. Additionally, in the case of the nMSSM, observable GW production is 
only possible (if at all) for small critical temperatures, what further reduces the prospects 
of chargino mediated baryogenesis due to Boltzmann suppressed chargino densities, even if 
they are not based on transport. 

Hence, the observation of GWs that can be attributed to the electroweak phase transition, 
could not only partially shadow a smoking gun signal from inflation, but also question (at 
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least the standard picture of) electroweak baryogenesis. 

It would be interesting to study baryogenesis with supersonic bubbles in more detail. 
After all, supersonic walls could appear in a larger part of the nMSSM parameter space. 
Finally, one should keep in mind that especially the turbulence contribution to the GW 
spectrum is still under debate. The available computations are based on non-relativistic 
models of turbulence. A numerical simulation of relativistic turbulence would be desirable 
to obtain more reliable predictions. Also the high frequency behavior of the contributions 
from bubble collisions is debatable. This was recently discussed using an analytic approach 



in Ref. [46| and will be subject of a forthcoming publication [47i] . 
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APPENDIX A: THE PHASE TRANSITION 

If the two minima of the effective potential are separated by a potential well, thermal 
and quantum fluctuations lead to small regions in space that acquire a finite Higgs vev close 
to the global minimum of the effective potential. If this region reaches a critical size, it 
is advantageous to expand this region, since the gain in energy by increasing the size of 
the bubble dominates over the increase in energy in the surface of the bubble. The phase 
transition proceeds then by nucleation and expansion of bubbles consisting of regions with 
non-vanishing Higgs vev. 



In the semi-classical theory of tunneling 



44 



45], the tunnel probability depends on the 



action of the so-called bounce solution. This configuration fulfills the boundary conditions 
denotes the symmetric minimum of the potential V{ip,T)) 

dM^) = 0, lim ^{g) = (Al) 

g— »oo 

and the Euclidean equation of motion 

^'^ + 2*£ = rM, (A2) 



dg^ Q dg 
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and 7 = 2(3) corresponds to tunneling at finite temperature (in vacuum). At finite temper- 
ature, the bubble nucleation rate is then given by 



where A is a constant of 0{1) and the Euclidean action is 

2 



(A3) 



5*3 = 47r / dg g 



(A4) 



To determine the bounce configuration, we use the method presented in Ref. 19[. This 
method is based on a two step procedure. First, the bounce solution is determined without 
damping, 7 = 0, which can be achieved relatively simple due to energy conservation. Sub- 
sequently, the parameter 7 is continuously increased to the desired value, 7 = 2, while a 
linearized version of the discretized equation of motion is iteratively solved. 

The temperature when the phase transition occurs depends besides the Euclidean action 
5*3 (T) on the cosmological parameters ji5|]. The expansion of the Universe is characterized 
by the Hubble parameter 

ff- = (A5) 



2^2 

where Mpi = 1.22 x 10^^ GeV denotes the Planck mass, and near the electroweak phase 
transition ~ 1/34. The probability that a bubble was nucleated inside the causal volume 
is given by 

(A6) 



-S3/T 



T 



The temperature of the beginning of the phase transition T„ is defined by the nucleation of 
the first bubble 

POO 

P\t=t. = dP = l. (A7) 

On the other hand, most gravitational radiation is emitted at the end of the phase transition 
when the bubbles collide. Neglecting the expansion of the Universe, the radius of a bubble 
at the temperature that nucleated at temperature T and expands with a velocity is 
given by [5| 



R{T^,T) = Vb- 



Thus, the fraction of the causal volume that is in the broken phase is 

AttH^ 



R\T^,T)dP. 



(A8) 



(A9) 
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We define the end of tlie pliase transition to be given by 

/It.=t, = 1- (AlO) 

Anotlier relevant quantity in tlie analysis of gravitational wave production is the typical size 
of the nucleated bubbles at the end of the phase transition, (R). Following Ref. {l^, we 
use the maximum of the bubble volume distribution, dV oc R^dP, which should be closely 
related to the typical scale for energy injection into the turbulent plasma. In this case, one 
obtains the condition (see Eqs. (I3l) and (1A8I) ) 

This definition has the virtue that it characterizes typical bubble sizes at the end of the 
phase transition in a reasonable way, even for very strong phase transitions, in contrast to 
P{Tf) that can become negative. The temperature Tr is by construction larger than Tf 
and typically smaller than T„. Notice that there is a certain arbitrariness in the choice 
of the length scale that enters into the analysis of the peak frequency of the turbulent 
spectrum. The authors of Ref. |l2| use the maximum of the bubble volume distribution in 



momentum space, given by {R) ~ 5vb/(3, while the authors of Ref. [13| prefer the size of 
the largest bubble (R) ~ Vb/p. We employ the maximum of the bubble volume distribution 
without transforming to momentum space what leads to Eq. (lAlip and compromises between 
these two choices. This arbitrariness is a major source of uncertainty in the GW spectrum 
contribution from turbulence. 

Finally, is the number of bubbles per causal volume at the end of the phase transition 

N = P\T=Tf ■ (A12) 

To obtain approximate results for these quantities, a Taylor expansion of the action can 
be used. Performing the integrals in Eqs. (1A7I) . ( 1A9I) and (1A12P leads to 



T 
T 



T=T„ 



141.4 - 4 log ( ^ ^ - log ( ) , (A13) 
^ ' lOOGeVy ^ V 100 ' ^ ' 



T=Tf 

3 



130.8-41og( — ^ ) -41og( ) +31ogt;6, (A14) 



(PiTfY) 
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